Iterative proportional fitting